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ABSTRACT 

We present pulse timing analysis of Rossi X-ray Timing Explorer (RXTE) 
observations of the anomalous X-ray pulsar IE 2259+586 from its 2002 out- 
burst to October 2010. Our objectives are to extend the work on the re- 
covery stage after the 2002 glitch, investigate the variations caused by the 
second glitch which occurred in 2007 and look for other unusual events, if 
any, that arise in the regular spin-down trend of the source. We find that 
the fractional change in the spin frequency derivative after the 2002 glitch 
is not stable since it decreased an order of magnitude, from —2.2 x 10~ 2 to 
— 1.278(3) x 10~ 3 , in about 2.5 years. From pulse timing analysis, we discover 
two small frequency shifts with fractional changes Avjv = 3.08(32) x 10~ 8 
and Avjv = —1.39(11) x 10~ 8 . While the first one of these shifts is not 
found to have a fractional frequency derivative change while second one has 
Av jv — —2.9(2) x 10~ 2 . We interpret these frequency changes as positive and 
negative microglitches similar to those seen in radio pulsars. 
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1 INTRODUCTION 



The X-ray pulsar (AXP) IE 2259+586 was discovered by lFahlman &: Gregory! (119811 ) within 



the supernova remnant G109. 1-1.0 with ~ 7 s p ulsations. Furt 



IE 2259+586 with 



1988 



Hanson et al. 



enm a , EXOSAT and Ging a (iKoyama et al 



1988 



Iwasawa et al. 



ler X - ray o b servations of 



19871 . 



1989 



Morini et al. 



19921 ) have indicated that the AXP is steadily 
spinning down with a rate of v ~ —1.0 x 10~ 14 Hz.s -1 , which is too slow to power the 



observed X-ray luminosity. Using ROSAT observations, iBaykal fc Swankj (119961 ) found a 
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marginal spin-up episode between 1991 and 1992. iMereghetti. Israel fc Stella! ( 119981 ) using 
the RXTE observations gave an upper limit on the se mi- axis to be . 03 It . sec suggesting a 
white dwarf companion smaller than 0.8 (see also 



observations, 



Baykal et al. 



(119981)). Using RXTE 
Kaspi. Chakrabarty fc Steinbergerl (119991 ) obtained, for the first time, a phase 



coherent timing solution for IE 2259+586 , indicating stability over a 2.6 yr time period. This 
stability strongly suggested the magnetar interpretat ion rather than the binary accretion 



hypothesis. Using ASCA and BBXRT observations, 



Corbet et al 



(119951 ) showed for the 



first time that the X-ray spectra can be explained by a blackbody plus power law model. 



INTEGRAL and RX' 



10 keV (IKuiper et al 



E ob servations have revealed the evidenc e of hard spectral tails above 



2006|) (see also Kaspi fc Bovdstunl (120101 ) ). 



In J une 2002, IE 2259+586 showe d an outburst consisting of over 80 short SGR-like 
bursts (jGavriil. Kaspi fc Woods! 120041 ) . During this outburst, flux level of the source in- 



accompanied with a g 



2003 



Woods et al. 



itch of fractional frequency 



creased by a factor of m ore than 20 
change ^ ~ 4.24 x 10" 6 (IKaspi et al. 

4r ~ 8.5 x 10~ 7 was ob served when the source was quiescent (iDib. Kaspi k, Gavriil 



2004). In 2007, a second g 



itch of 



Woods et al. 



20081 ) 



(120041 ) modelled XMM and RXTE spectra of IE 2259+586 before and 
after the 2002 outburst with the blackbody plus power law model modified by interstellar 
absorption. The blackbody temperature increased rapidly after the outburst in contrast 
to the photon index which decayed about half of its pre-outburst value. Both parameters 
quickly recovered t o within 25% of their their pre-burst levels within the first ~l-3 days 
( Woods et al. 2004^ . In the first few hours of the outburst, XMM and RXTE observations 
revealed that pulsed fraction decreased to ~ 15% then, like the spectral parameters, it 
quickly recovered to its pre-outburst state within 6 days. There were also some variations 
in the phase dependence of the energy spectrum: While the photon index showed significant 
variability one week before the outburst, this variation vanished three days after the outburst. 
XMM-Newton observations of the source showed that X-ray flux decayed with a power law 
index ~ —0.69 after 2002 ou tburst unti l 2005 . A strong correlation between X-ray flux and 
hardness ratio is observed by IZhu et al.l ( 120081 ) . Pulse profile of the source changed in such a 
way that much of the power moved to the fundamental harmonic and the second harmonic 
regained i ts power within 1 w eek after which a very slow recovery within seve ral weeks was 



observ ed ( iWoods et al.ll2004j ). Using the cumulative properties of the outburst. IWoods et al. 



( 120041 ) concluded that the pulsar experienced something which was not sudden at all and 
including two components, one of them was on the surface such as a series of fractures 



and the other distributed over a wider region like smoother plastic change. Hence, both the 

superfruid interior and the magnetosphere were af fected. 

The Spitzer Space Telescope observations of Kaplan et al.l (120091 ) at near- and mid- 



infrared regions of the electromagnetic spectrum w as interpreted as a p ossible indication 



of the presence of a passive X-ray-heated dust disc. 



Kaplan et al.l (120091 ) also showed that 



the IR data of this source can also be modelled by a power law spectrum which could be 



explained by magnetospheric emission. 



Tian. Leahy fc Lil ( 120101 ) estimated the distance to 



the AXP (and also to the SNR G109. 1-1.0) as 4.0 ± 0.8 kpc using the 21-cm H i-line and 
CO-line spectra of the SNR G109. 1-1.0, H n region Sh 152, and the adjacent molecular cloud 
complex. 

In this paper, we construct phase coherent pulse arrival times between 2002 and 2010. 
Using these arrival times we refine the timing solution after the 2002 glitch and look for 
frequency and derivative changes in these time intervals. 



2 DATA & ANALYSIS 
2.1 Data 

The dataset used in this work consists of Rossi X-Ray Timing Explorer (RXTE) proportional 
counter array (PCA) observations of the AXP IE 2259+586 covering the time between March 
2000 (MJD 51613) and October 2010 (MJD 55483). The PCA operates in the energy range 
2-60 keV using an array of five collimated xenon/methane multi-anode proportional counter 
units (PCUs). T he instrument has a total effective area of ~ 6500 cm 2 and a field of view 



of ~ 1° FWHM (IJahoda et al 



19961 ) 



A total of 479 observations were used for the analysis presented in this paper. The 
durations of observations vary from 10.6 ks to 0.2 ks, most of them being greater than 1 ks. 
During the analysed RXTE-PCA observations, the number of active PCUs varied between 
1 and 4. Due to timing concerns, we used all the available layers of all PCUs in our analysis. 



2.2 Pulse-Timing Analysis 

The GoodXenonWithPropane or GoodXenon mode and event mode data were energy se- 
lected (2-10 keV) for all xenon layers, and were binned with 125 ms time resolution. The 
time values in the light curves were then corrected to the solar system barycentre. We used 
the standard analysis tools for RXTE-PCA data included in FTOOLS package to obtain the 



to avoid sud- 



Woods et al. 



4 B. Igdem et al. 

light curves and to merge them into a single light curve covering all the interval mentioned 
above. In the extraction of lightcurves we removed all outburst events in order 
den p ulse profile changes. The light curves were folded at frequencies given by 
( 12004] ) . The time series was split into intervals of approximately same duration, and each 
segment was folded with a quadratic ephemeris with the same frequency and frequency 
derivative, so that we obtained a pulse profile for each time interval, which is made up of 20 
ph ase bins. Then, we switched to the harmonic representation of pulse profiles as introduced 



bv lBovntonl fjl985l ) 



10 



F o + 22 F k cos k((f) - <j) k ), 



(1) 



fc=i 



The template pulse is extracted from a longer time interval at which glitch recovery has 
already taken place 



10 

g{cj>) = ^GfcCosA;^-^). (2) 

k=l 

(3) 

Then, by cross-correlating the pulse profiles with the template pulse, we obtained the pulse 
arrival times, A<p. In order to see the effect of pulse profile changes in pulse timing, we also 
performed pulse timing using first and second harmonics independently. Our pulse arrival 
times were consistent within la level with 10 harmonic expansion of pulse profiles for 20 
phase bins. This suggests that even if the pulse profile changes, pulse phases are not effected. 
As a final check, we gave weight to all pulse harmonics using the variances of harmonics 
from all pulses. We applied the cross correlation for the weighted pulses (or filtered pulses), 
and again we obtained the pulse arrival times as those of unweighted pulses in la level 
confidence. We concluded that pulse shape changes do not have any effect on pulse phases 
and pulse timing analysis. 

In top panel of Fig. [fl it is seen that the slope of the pulse arrival times changes after the 
first glitch. In the estimation of pulse arrival times, the spin down rate of the source was kept 
fixed at the value before the glitch, which is v = —9.92 x 10~ 15 Hz.s -1 . By modelling the 
pulse arrival times lying after the first glitch by a first order polynomial, with A0 = fo + Sist, 
the new pulse arrival times were obtained for the corrected pulse frequency, u + 8u. In the 
second panel of Fig. [TJ it is clearly seen that the slope of the pulse arrival times is shifted 



again after MJD 54300. The shift in the pulse frequency is found again using the method 
explained above and new pulse arrival times are obtained after MJD 54300 by correcting the 
pulse frequency again. The third panel of Fig. (TJ where we present pulse arrival times for 
the corrected frequencies after the first and second glitches, shows that there are additional 
breaks in pulse arrival time series around MJD 53750 and 54880 (see vertical dash lines 

quoted B and E) 

The first glitch is at MJD 52443.13 and was reported by lKaspi et al.l (120031 ) to be accom- 
panied by an outburst of many individual bursts. About 1360 days after the first glitch, at 
around MJD 53750, we resolved a small fractional frequency shift of = 3.08(32) x 10~ 8 . 
Phase connected pulse arrival times after MJD 54300 indicates another glitch with a frac- 
tional frequen cy change o f ^ ~ 8.20(2) x 10~ 7 . This second glitch is consistent with 1 a 



level given by 



DibetaL 



( 120081 ) in their review paper of AXP glitches. Moreover, around 
MJD 54880, we have observed another fractional frequency shift of = —1.39(11) x 10~ 8 
with = —0.029(2). All the glitch episodes in this work are presented in Fig. [1] and timing 
solution are given in Tables [1] and |2j 



For the first glitch, as described in 



Woods et al 



(120041 1 . the standard post-glitch relax- 



ation models, which consist of single or multiple exponential relaxation terms (Alpar et al. 
1984b), did not provide reasonable fits. An alternative model is the irregular model includ- 
i ng an exponential rise term, which had been developed for two gli tches of Crab pulsar 
(ILyne. Pritchard Graham-Smith! Il993l ; IWong. Backer fc Lynd l200ll ). In this model, the 



increase in the slowdown rate is interpreted as the cumulative of successive glitches; part of 
the increase is provided by a step and part is in the form of an exponential rise 



v = v + i>o(t - to) + Az/ + Au g (l - e~ (t - ts)/Ts ) 

-Au d (l - e -(*-*s)/^) + Aut, (4) 

where the first two terms represent the frequency evolution before the glitch, t g is the glitch 
epoch, Av is the frequency jump with the glitch, Av g and Av d are the growth and decay 
amplitudes, respectively, r g and t& are the growth and decay time- scales respectively, and Av 
is the jump in the spin down rate. We used the corresponding phase evolution equation to fit 
our pulse arrival times in the period MJD 52390-53750 (The correction to the pulse frequency 
after the first glitch mentioned above is 5v = Au+ Av g — Av^). In the fitting process, we used 
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Table 1. Timing Solution for IE 2259+586 Before MJD 53750 



Parameter 


Value 1 


Value 2 


Spin frequency,^ (Hz) 


0.14328703257(21) 


0.14328703257(21) 


Spin frequency derivative, u (Hz.s -1 ) 


-9.920(6) x 10- 15 


-9.920(6) x 10- 15 


Epoch (MJD) 


52400 


52400 


Au (Hz) 


5.25(12) x 10~ 7 


6.70(1.14) x 10~ 7 


Au B (Hz) 


> 8.7 x 10- 7 


8.29(78) x 10~ 7 


T a (days) 


14.1(7) 


14.1(1.2) 


Av d (Hz) 


Av g +(~5x 10~ 9 ) 


1.06(8) x 10~ 6 


r d (days) 


15.9(6) 


15.9(1.1) 


Ai> (Hz.s" 1 ) 


2.18(25) x 10~ 16 


1.268(3) x 10" 17 


tg (MJD) 


52443.13(9) 


52443.13(10) 


rms timing residual (ms) 


44.9 


85.1 


Start observing epoch (MJD) 


51613 


51613 


End observing epoch (MJD) 


52900 


53750 



Numbers in parentheses indicate la uncertainti es in t he least significant digits quoted. 1 Taken from table 4 of lWoods et al.l 

d2004h . 2 From our analysis. 

the x 2 statistic and leave all parameter s as free parameter. For normalizations of exponentials 



Woods et al. 



(120041 ) . 



we used the initial ranges as given by 

In Table (U we present timing solution of the first glitch for the extended time coverage 
from MJD 51613 to MJD 53750. Our analysis indicates that 4^ has shown a variation 
from —2.2 x 10~ 2 to —1.278(3) x 10 -3 , i.e. there is an order of magnitude reduction in the 
fra ction in a ~ 2 . 5-yea r time. We also fitted the pulse arrival times for a time interval given 



by 



Woods et al. 



( 120041 ) and found that the timing solution is consistent with Woods et al. 
(2004) in 1 a level. 

The small frequency shifts of ~ 10~ 8 with both signs at MJD 53750 and 54880 are 
at the order of microglitches seen in radio pulsars. In the first frequency shift, ^ is not 
significant, however in the second one this value is found to be —0.029(2), i.e. the fractional 
change in the pulse frequency derivative is negative. This is not very unusual when compared 
to the microglitches in ra dio pulsars, since there is no preferred sign for both of the fractional 
jumps Avjv and Az>/z> (jChukwude fc Uramall2010l ). The timing parameters of these small 
frequency shifts are given in Table |2j 



3 DISCUSSION 

In this work we presented extended pulse timing analysis of IE 2259+586 from 2002 to 
2010. In phase coherent pulse timing analysis, we found that after 2002 glitch reduced 
an order of magnitude in a ~ 2. 5- year time. We also discover two micro glitches Avjv = 
3.08(32) x 10~ 8 and Avjv = —1.39(11) x 10~ 8 . In the first frequency shift, we do not resolve 
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Table 2. Timing Solution for IE 2259+586 After MJD 53750 



Parameter 


Microglitch 1 


Glitch 2 


Microglitch 2 


Spin frequency,!/, (Hz) 


0.143286381(13) 


0.143286138(14) 


0.143285760(14) 


Spin frequency derivative, v (Hz.s -1 ) 


-9.920(6) x 10" 15 


-9.920(6) x lO" 16 


-9.920(6) x 10" 15 


Epoch (MJD) 


53750 


54040 


54880 


Au (Hz) 


4.42(46) x 10~ 9 


> 1.6 x 10~ 6 


-2.00(15) x 10- 9 


Ai> (Hz.s- 1 ) 






2.91(22) x 10~ 16 


tg (MJD) 


~ 53750 


~ 54040 


~ 54880 


rms timing residual (ms) 


76.8 


66.8 


106.7 


Start observing epoch (MJD) 


53700 


53900 


54800 


End observing epoch (MJD) 


54040 


54900 


55570 



Numbers in parentheses give lcr uncertainties in the least significant digits quoted. 

frequency derivative. In the latter one, fractional change in the spin down rate is found to 
be Au/u = -2.9(2) x 1(T 2 . 

Microglitches make up a class of small amplitude resolvable jumps tha t can be observed 
in both, or either of pulsar rotation 



1985 



Chukwude 



2002 



Cordes et al 



requency or its first time deriv ative ( jCordes fc Downs 



1988 



D'Alessandro et al 



19951 ). The orders of magni- 
and 10" 5 < |Az>/z>| < 10~ 2 



tude for v and v are in the range 10 < \Avjv\ < 10 
respectively. So far 319 m i cro g li tches are observed from 46 radio pulsars ( jCordes fc Downs 



1985 



D'Alessandro et al 



1995 



Chukwude fc Urama 



20101 ) . Microglitches observed in IE 



2259+586 are of the order of the upper limits of radio pulsar microglitches in magni- 
tude . Unlike the large glitches, microglitches can have both negative and positive signs 
in their fractional pulse frequencies and first time derivatives. Henc e,it is difficult to under 



stand microg 



Cordes et al 



itche s in terms of theories developed for large glitches ( jCordes fc Downs 



1985 



19881 ). The cause of micro glitches is unclear but can b e explained as involvin g 



neutron star interior and/or magnetospheric torque fluctuations (IRuderman et al.lll998l ). 
he magnetar crust ca n be cracked by the variations of the magnetic field configuration 



( IRuderman et al.lll998l ). In this case, both signs of microglitches may be observed. 



A nti-glitch (negative g litch) behaviour have been seen in two other magnetars ( jThompson et al 



2000 



Gavriil et al 



201 if) . Sudden spin down trend of SGR 1900+14 is interpreted by 
Thompson et al.l (120001 ) as an anti-glitch with a fractional pulse frequency change of Avjv ~ 
10~ 4 . However, the scale of this event is four orders of magni tude larger than that of the 
negative microglitch of IE 2259+586. For AXP 4U 0142+61, Gavriil. Dib fc Kaspil ( 2011 ) 
showed that the total fractional frequency jump at the time of glitch is Avjv ~ 1.7 x 10 -6 ; 
however, the decay of large glitch event left a net negative slow-down portion of Avjv ~ 
—8 x 10~ 8 after 17 days of recovery. In addition, for rotation powered pulsar PSR J1846- 



8 B. Igdem et al. 

0258. lLivingstone. Kaspi fc Gavriill (120101 ) reported a large spin up glitch which experienced 



a large recovery and ended up with a net decrease of pulse frequency similar to the recovery 
of AXP 4U 0142+61 and event occurred in SGR 1900+14. So far IE 2259+ 586 is the uniqu e 
source showing microglitch/ anti-glitch events among the AXPs and SGRs (IDib et al.ll2008l ). 



In magnetars, glitches may be triggered by str ong internal magnetic 



fields as the crust 



is deformed, either plastically or cracked violently ( iThompson fc Duncan 



19961 1 . There may 



be some parts of superfluid rotating more slowly than the crust. Then, by a glitch the 
angular momentum is transferred from the crust to the superfluid, opposing the case in 
ra dio pulsar glitches , and t he crust spins down after the glitch. This explanation is proposed 



by 



Thompson et al.l (120001 ) to explain the net spin down rate event in SGR 1900+14. Hence, 



small scale magnetar glitches can explain the glitch/anti-glitch trend of micro glitches like 
those seen in IE 2259+586. 

The 2002 glitch is well-studied by 



Woods et al. 



( 120041 ) until the post glitch recovery 



mid of 2003. We extended post glitch recovery until the first micro glitch event seen in 
2006. In the pulse timing analysis, we confirmed that the fractional frequency change of the 



glitch has an amplitude of Au/v 



Table 1) determined by 



Woods et al. 



. 68(80 ) x 10 6 , and other post glitch parameters (see 
(120041 ) from their 500-day time span of observations. 



The fractional change in the pulse frequency derivative Az>/z> = —1.278(3) x 10 3 , from 



Woods et al. 



our 13 60-day pulse arrival times is found to be 17 times smaller than that of 
( 120041 ). which is deduced from a shorter time interval. We should note that our analysis in 
their timing interval yielded a similar value in 1 — a limit. (This 1 — a consistency in Az>/z> 



does not disru pt even if we 
values given by lWoods et al. 



reeze the normalizations of the exponential components to the 
( 120041 ) ). This indicates that f ractional change of pulse frequency 



derivative is variable. As discussed in 



Woods et al. 



( 120041 ) the negative sign of the fractional 



change of derivative of pulse frequencies is quite unusual for large radio pulsar glitches. For 
the large glitches fractional frequency change is Au/u > 10~ 7 and fractional pulse fre quency 
derivative has the range Az>/z> ~ 10~ 2 - 10^ 3 (see also lAlpar fc Bavkall (11994 120061 )). 

In the standard large glitches of radio pulsars, the sudden shifts in the pulse frequency and 
the spin down rate of a radio pulsar have been considered to be due t o the coupling mecha- 
nism between the crust and the crust superfluid (Alpar et al. 1984a,b); (lAlpar. Cheng fc Pines 



19891 ). The angular momentum transfer between the crust superfluid and the crust is pro- 



vided by the quantized vortices. There are two types of vortex lines separated according to 
their dynamical behaviour in the superfluid. One type is responsible for the continuous vor- 



tex flow since they are not pinned to the inner crust of the neutron star but they unpin and 
repin at some rate. The other type of vortices are pinned to the crustal nuclei and are not 
allowed to move radially outwards up to the maximum tolerance of the pinning force. Hence, 
a lag appears between the spin-down rates of the crust and the superfluid, v = v s — v c > 0. 
When the vortices cannot stand this lag any more, they suddenly unpin from the crustal 
nuclei and some angular momentum is transferred to the crust. This sudden unpinning of 
the superfluid neutron vortices from the crustal nuclei can result in glitches of fractional 
pulse frequency changes of A1//1/ ~ 10~ 7 to 10~ 6 together with increases in the magnitudes 
of spin down rate. The latter type of vortices do not have any contribution to the usual spin 
down trend of the pulsar other than the glitches. Hence, the relaxation of the spin down rate 
after the glitch is accomplished by the continuous vortex flow, and the persistent change 
in the spin down rate, which is observed sometimes, is interpreted as the increase in the 
amount of the pinned vortices that cannot move towards the outer crust after the glitch. If 
the superfluid part supplying the vortex flow has moment of inertia I a and the trapped part 
has the moment of ine rtia Jr, then the equations governing t he dynamics of the glitches are 
determined as follows (lAlpar. Kiziloglu fc van Paradijslll995l ) 

'I_A I_B 

21 I 



j 

V 



(5) 
(6) 



Au 

V 

Az> I_a 

V I 

where I is the total moment of inertia of the pulsar and Su is the change of the rotation 
frequency of the superfluid interior. Thus, in the vortex creep theory the fractional change 
in the spin down rate is determined by the fraction I a/ I which should have a positive value 
in the presence of a constant external torque applied on the pulsar. Unless the vortices prop- 
agate inwards to the rotation axis of the star, the negative value of 4^ can be e xplained by 



exter nal torques. This could be due to the secular change of magnetic moment ( iRuderman 



199ll ) or change in the electromagneti c torque by a reorgan ization of magnetic moment and 



angular momentum axis of the star (ILink &: Epstein! 119971 ). However, fractional change of 
the pulse frequency derivative by a factor of 17 in a time scale of only 860 days due to the 
magnetic moment change is quite fast. Another ex planation could be rela xation of twisted 



magnetic fields which creates stresses on the crust (jThompson et al.l 120021 1 . Persistent seis- 



mic activity may modulate the spin down to rque and crustal shear waves drives particles 



wind and causes angul ar momentum changes (jThompson fc Blaes 



1998 



Harding et al. 



1999; 



Thompson et al.l 120001 ) . It is interesting to note that all glitches of IE 2259+586 have oc- 
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curred almost at a constant spin down rate which is i> = —9.920 x 10~ 15 Hz ~ l . This suggests 
that glitches happen to be when post glitch frequency derivatives are completely relaxed. 
Therefore, the excitation of precession during the glitch, due to a stable reorientation of 
magnetic axis relative the rotat ion axis of the star is considered unlikely since long term 



torque on the star is not stable (ILink fc 



As a final remark, after 2002 glitch, 



.ostein 



1997 



Kaspi et al.l (120031 ) report Az>/z> ~ 1.1 for a time 



span of a couple of weeks. They suggest that the increase of spin down rate is due to the 
decoupling of crust and core superfluid where they re couple in 2 weeks time-scale which is 



extre mely long compared to the theoretical estimates (lAfpar fc Sauls 



1988 



Sidery fc Alpar 



20091 ). This is also clearly resolved in our pulse timing analysis for a particular short time 



span; however, as we are interested in new events in the long term recovery and unusual 
events in pulse timing, this aspect of 2002 glitch is kept out of scope of this work. 
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Figure 1: Phase offset series for AXP IE 2259+586. Panels are described from top to bottom. 
Top panel: Phase offsets extracted using the spin-down model of the period before 2002- 
glitch. Second panel: Phase offsets using the correction of the pulse frequency after MJD 
53443.13. Third panel: Arrival times obtained by an additional correction after MJD 54300. 
The solid curve is the phase offset evolution of the models presented in Tables [1] and |2j 
Bottom panel: Residuals, after subtracting the best-fitting models given in Tables Q] and |2j 
All panels: Glitch epochs are indicated with dashed vertical lines. A: 2002 glitch corresponds 
to MJD 52443.13, B: small microglitch at MJD ~53750, C-D: no phase-connection, D: second 



